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I.  INTRODUCTION 


Optimization  and  synthesis  methods  for  antenna  arrays  have  been 
studied  by  many  researchers;  the  exi stance  alone  of  a long  series  of 
papers  on  these  subjects  is  enough  to  emphasize  the  importance  of  these 
areas.  Specifically  Dolph  [1]  and  Riblet  [2]  have  published  a series 
of  very  interesting  papers.  These  papers,  by  using  Chebyshev  polynomials, 
give  a synthesis  for  uniform  linear  arrays  which  offers  the  minimum  beam- 
width  for  a prescribed  sidelobe  level.  Unz[3-6],  Harrington[7] , King, 
et  al.  [8],  Jacabsen  and  Madsey  [9]  and  many  others  gave  answers  to 
very  interesting  problems.  On  the  array  optimization  Uzkov  [10]  first 
found  by  using  linear  transformations,  the  maximum  directivity  of 
linear  arrays.  Uzkov's  theoretical  work  was  extended  by  Bloch,  et  al . 
[11],  Uzsok  and  Solymar  [12]  and  Stearns  [13].  Tai  [14]  considered  the 
problem  of  achieving  maximum  directivity  in  uniform  linear  arrays  of 
short  dipoles  and  gave  many  useful  graphs  on  the  subject.  Cheng  [15] 
studied  the  determination  of  directivity  in  more  complicated  arrays. 

Lo,  et  al . [16]  obtained  the  optimum  SNR  with  constraint  on  the  Q-factor. 
Sungiri  and  Butler  [17]  have  used  the  eigenvalue  method  to  find  the 
solution  for  the  maximum  directivity  with  constraints  on  the  resulting 
sidelobes.  Matrix  methods  were  applied  by  Strait  and  Kuo  [18],  Sarkar 
and  Strait  [19]  and  Sahalos  [20]  for  constrained  optimization  of  various 
performance  indices  of  arrays  with  straight  parallel  or  nonparallel 
thin  wire  antennas. 

The  orthogonal  method  for  arrays  consisting  of  arbitrarily  oriented 
short  dipoles  was  applied  by  Sahalos  [21-22].  The  orthogonal  method 
one  decade  ago  was  used  first  by  Unz  [23]  and  recently  by  Sahalos  [24-26] 
in  many  antenna  array  applications. 

In  the  present  work  an  effort  is  made  to  give  some  useful  formulas 
applicable  to  synthesis  and  optimization  problems  when  the  arrays  con- 
sist of  nonparallel  wire  antennas. 


II.  MOMENT  METHOD  FORMULATION 

Let  a wire  structure  (Figure  1)  be  composed  of  a number  of  straight, 
segments.  By  defining  a right  handed  orthogonal  coordinate  system  (n,,?,'') 
at  each  point  of  the  wire's  cylindrical  surface,  we  can  find  [27]  that: 

- . (1) 

This  expression  expresses  the  reaction  integral  equation  developed  by 
J.  Richmond  and  is  true  for  electrically  thin  wires. 
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By  expanding  the  current  on  the  wire  in  a finite  series  of  the 

form 


1(0  = I I/„(0  (2) 

n 

we  can  have  from  Equation  (1)  a system  of  simultaneous  linear  algebraic 
equations 

[Z](I)  = (V)  (3) 


where  the  elements  of  the  matrix  [Z]  are: 

The  quantities  in  Equations  (1)  and  (4)  are  given  by: 


‘-m 


ni 

(Ji  • E^  - • i^f)  dv 

(5) 

i_ 

r^T  ^ ^ 

« • E*"  d(t 

(6) 

Ztt 

0 

1 

2t- 

f2T.  * 

<(>  • d4. 

(7) 

where  (J^,M^).ar^  the  impressed  currents  by  which  the  wire  structure  is 
excited  and  (E'",H"’)  is  the  field  of  an  electric  test  source  located  in 
the  interior  of  the  wire  surface  and  radiating  in  free  space. 


By  the  help  of  Equation  (3)  \/o  can  find  a relation  between  the 
currents  and  the  corresponding  voltages  of  the  fonn; 


(I)  = CY](V) 


(8) 


where  in  the  column  (V)  the  elements  are  nonzero  only  at  the  feed  points. 
The  far  zone  field  transmitted  by  an  antenna  can  be  expressed  as: 


E((|>,0)  = £*’(<(>,6)  + £'^(41,0) 


(9) 
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If  the  wire  structure  is  composed  of  N1  segments  and  N(N  < N1 ) main 
ports,  the  expression  of  Equation  (9)  will  become: 

, N r N1  1 

E{>i.e)  = I V J I [E®{  Ke)  + En^.o)]^  . (10) 

i=l  n=l  " 

^ - V 

The  vectors  E{^  and  Ej^  are  the  electric  field  components  of  the  n^h 
segment.  If  the  center  of  the  segment  has  coordinates  and 

direction  angles  then: 


E^(.f,9) 


[cose  sine^  cos((i)-ij)^ ) - sine  cose^]  exp[jq^-]o 'j 
sine^  sin(!!>’-(t.)  exp[jq.]4,  f 


where: 


q^  = ^r^[sine  sine^-  cos(({'-4i^. ) + cose  coso^]  . (12) 

III.  ORTHOGONAL  METHOD  FORMULATION 

Equation  (10)  shows  that  the  electcic  fiel,d  is  a vector  of  an  Nl- 
dimensional  vector  space  with  a basis  {Oil  = ^E^'q  + ?|}.  As  we  can  see 
N vectors  of  this  space  define  the  E((t>,ej.  By  a method  analogous  to 
that  of  orthogonality  we  can  express  the  field  as  a situation  of  N ortho- 
gonal vectors  of  the  same  space.  Let 


^‘i  = I Y„i(?S  . eJ) 

n=l 


♦j  ” 1 ^ ■ C'') 

m=l 

The  inner  product  of  these  two  vectors  is: 

r N1  N1  rn  ,2tt 

‘**’*J’  ° Ul  m.^1  ''"'■''"J  Jo  lo 


We  suppose  that 
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The  expression  of  can  be  found  as  a function  of  the  array  geometry 
and  is  given  in  Appendix  I.  The  electric  field  can  now  be  expressed 
with  the  help  of  'ifj  as: 

N 

E((f.o)  = Lj  ii)j((l),9)  (24) 

i = l 

where  the  L^-  are  given  by 

(25) 


and  the  corresponding  voltages  by  the  help  of  Equation  (20)  will  be 


Vi  = .1. 


J = i 


(j) 


(26) 


As  result  from  the  above  discussion  the  synthesis  of  an  array  will 
follow  these  steps: 

i)  Definition  of  the  five  coordinates  of  each  segment  (three 
positional  and  two  directional); 

ii)  Calculation  of  the  constants  C; 

iii)  Evaluation  of  the  field  E(ij).o); 

iv)  Calculation  of  ; and 

v)  Computation  of  the  feed  voltages  V^- . 
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IV.  OPTIMIZATION  PROBLEM 


The  optimization  problem  is  a procedure  of  maximization  of  an 
index  of  the  array  with  or  without  constraints.  Some  of  the  usual 
performance  indices  which  are  involved  follow. 


1.  Power  Gain 


G = 4^  radiation  intensity  for  specified  direction  . (27) 

power  input  to  the  array 


In  relation  to  the  feed  voltages  and  the  input  admittances.  Equation 
(27)  is  expressed  as 


N , 1 2 

y V-  ()).(  i ) 

1 ' *^0  0\ 


N N 

15  y 1 ^ '^ii^ 

i=l  j=l  ' J ’J 


(28) 


where  Y^j  is  the  admittance  between  the  main  ports  i,j. 


2.  Directive  Gain  and  Directivity 
The  directive  gain  is  defined  by 


D = 4tt 


radiation  intensity  for  specified  direction 


radiation  power 

It  can  be  shown  that  the  radiation  power  is  an  expression  of  the  form: 


N N 

p,.  = y y v.v*(y';‘.  + y**?  - y 

i=l  J 


i=i  jri  1 3’  IJ  13  13  13' 


(30) 


where  complete  details  for  the  Y^j  are  available  in  Reference  [28], 


The  expression  of  directive  gain  becomes 
,I  V, 


i = l 


(31) 


and  the  maximum  value  of  D(('QnQ)  is  the  directivity. 
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3 . Efficiency  Indices 


The  efficiency  indices  are  very  important  parameters  describing 
the  performance  of  an  array.  Many  indices  can  be  defined  by  several 
ways.  Two  of  these  are 


2 _ radiation  intensity  in  the  direction  of  max,  radiation 
sum  of  the  excitation  voltages  magnitudes  squared 


squared 


5 _ radiation  intensity  in  the  direction  of  max,  radiation  (33) 

1 “ sum  of  the  feed  port  current  magnitudes  squared 


Equations  (32)  and  (33)  can  be  written 

N 

J,  ''i 

C _ I 1 =1 


= 


i^l 


N 

Vi 

1=1 


N 

I 

i = l 


N 


I ViY.  . 
j=l  ' ’J 


K = and  t (sensitivity  factors) 


(34) 


(35) 


Some  other  factors  as  sensitivity  and  Q-factors  can  be  defined.  These 
are  related  to  the  above  indices  and  given  by 


(36) 


c c 

Q = 5 and  Qi  = 5^  (Q-factors), 


(37) 
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A . Gain-Directive  Gain  and  Efficiency 
Indices  Maximization 

Until  now  we  have  seen  that  all  the  above  indices  can  be  written 
in  the  following  formula 


.1  ^'i('^o’^o) 


In('^o*«o)  - Sr-li ^ 

y V v-v^s-- 

i^i  j=i  ’ ^ 

where  i is  expressed  correspondingly  in  Equations  (28),  (31),  (34) 
and  (35). 

By  the  orthogonal  method  in  the  denominator  of  Equation  (38)  we 
can  find  an  expression 


N M N 

I I ViV^S.^  = I L . (39) 

i=l  j=l  1 i=l  1 T 

The  L.  are  related  with  V.,-  by  Equation  (26).  The  necessary  factors 
c('^Ure  given  from  Equation  (23)  by  substituting  the  K.jj  by  S^j.  Using 
the  expression 

Equation  (38)  yields 

IN,  2 

I I L.  ■i’.(T  ) 

I ^ . (41) 

‘-iLT 

1=1 

we  have 
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I ( '(>  , G ) — 
n' ^o’  o' 


N 

2 

2 

.1  Li  , 

+ 

.1  h » 

i = l 

J2J 

N 

y L-Lt 
i = l ^ ^ 


(42) 


where  = ip.  • o and  By  the  Schwartz  inequality  we  can 

take  ’ ^ ’ 


" L ■,<> 

‘-i  ‘i 

1 = 1 


/IN  \ / N \2 

(J, 


N 


e .*0 


N 

y L.L* 
i = l ’ ' 


N 

y L.L* 

i=l  ^ ^ 


4,  ‘'i** 


(43) 


and 


N 

2 /N  \/^ 

.1  h *1 

1 = 1 

(J, 

N 

- N 

.1 

1 = 1 

I L.L^ 

i = l ^ ’ 

From  Equations  (43)  and  (44)  we  can  see  that 
^ .0  *G 


i = l 


(44) 


n '^vi  ' 

1 = 1 


By  the  same  procedure  as  in  Reference  [22]  we  can  show  that 
Li  = k *0*0  + A 

where  the  parameters  k and  A are  related  by 


(45) 


(46) 
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l-K'^ 


N 

I 

1=1  ’ ^ 


= N 


1 = 1 


(47) 


Equation  (47)  gives  at  the  same  time  through  K and  A the  polarization 
sought. 


B.  Optimization  Subject  to  Constraint 
on  Electric  Field  Nulls  and 
Sidelobe  Levels 


Suppose  that  we  want  to  determine  the  feed  voltages  that  will  pro- 
vide electric  field  nulls  and  sidelobe  levels  in  K directions  for  the 
E^  component  and  in  \ directions  for  the  E*  component.  The  above  con- 
straint give  K+a  relations  between  the  feed  voltages.  If  the  field  in 
the  K direction  is  related  with  the  Eq  in  the  direction  of  maximum  by 
the  form: 


m 1 

Kl' 

E^l 

K2 

■ 

X 

• 

• 

• 

1 

(48) 


and  correspondingly  for 


r ~i 

! Et 

Al] 

E| 

A^i 

. 1 

1 

• 

1 

. 1 

• 

. 

Av 

A_ 

. >d 

E<t> 

0 


the  eJ  field 


= [A]EJ 


(49) 
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then  by  the  help  of  the  Equation  10)  we  can  find  a relation  between  the 
K+'.  feed  voltages  and  the  rest  N-K-a,  More  details  for  this  relation 
can  be  found  in  Reference  [20].  Finally  we  will  have  a modified  expression 
of  In('l’o’*^o)  which  will  be  of  the  form 


N-K-A 

V.  2.(,„.6^) 


N-K-A  N-K-A 

I I V^V*g 
i=l  j=l  ^ J ’J 


(50) 


Equation  (50)  is  similar  to  Equation  (38)  and  the  maximization  of  the 
index  will  be  done  by  the  same  method  as  before. 

The  g^  • are  the  matrix  elements  [Q]  of  Equation  (58)  and  the 
are  the  elements  of  the  vector 


[Z]  = [r^lo  + [E2]* 
of  Equation  (56)  in  Reference  [20]. 


The  above  procedure  gives  the  maximum  index  1^  when  K+A  values  of 
the  electric  field  are  known. 

If  we  want  some  of  these  to  be  the  si  delobe  levels  then  we  can 
use  an  iterative  procedure.  Thus  we  find  first  the  feed  voltages  by 
the  above  procedure  and  the  level  of  the  electric  field  sidelobes. 
These  are  compared  with  the  desired  levels.  The  new  directions  of 
sidelobes  are  used  to  find  new  feed  voltage,  etc.  The  procedure  is 
continued  until  all  the  sidelobes  take  on  the  desired  values. 


V.  EXAMPLES 

A.  Synthesis 

In  all  examples  we  assume  that  the  wires  have  radius  0.005a,  are 
centerfed  and  their  length  is  a/2.  For  the  case  of  synthesis  a double 
integration  is  applied.  In  Appendix  II  we  show  a method  of  numerical 
integration  with  good  accuracy.  This  method  is  extended  in  Reference 
[30]  and  uses  the  Chebyshev  polynomials. 

1.  Let  us  suppose  at  first  that  we  have  two  dipoles  normal  to 
each  other  with  their  centers  a distance  about  equal  to  the  diameter 
of  the  wire.  We  wish  to  have  an  electric  field  of  the  form 

f((j>,9)  = eJ'*’  (<. 


) 

i 
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By  the  orthogonal  method  we  found  the  feed  voltages 


= 1.  and  V2  = -j 

The  above  values  give  the  well  known  omnidirectional  dipole. 

2.  Suppose  a directive  field  of  the  form: 

E(iii,0)  = exp(-a-|Sin^e  cos^ii)  cos(2bisin0  cos4’)0  + 

00 

+ exp(-a2Sin^e  cos^i}')  cos(2b2Sine  cos4>)<t> 
has  the  following  properties: 

The  E®  is  max.  when  sinecos<j)=0,  zeroing  of  the  E®  mainlobe  is  when 
sinecoS(()=Tr/4bi  and  the  beamwidth  of  E®  is  2sin0^cos(}>i  when  /2/2=exp(-af 
sin^o^costtc]  )cos(2b')sino^cos((i-| ) . In  the  same  way  we  can  find  the  Ei*. 

An  array  with  N=7  dipoles  in  equal  distance  x/2  and  with  directions 
41  =0,  0'=ii/4  as  shown  in  Figure  2,  for  a'|=a2=2.683  and  bi=b2=1.428  gives 
the  following  feed  voltages  \/i=4. 237/11°,  \/2=V3=3. 735/5° . V/|=V5=1 .76/8°. 3 
and  V5=V7  = .27/aU.. 

If  Figures  3a, b we  can  see  the  normalized  E®  and  E'**  field  as  a 
function  of  p=sin0  for  (j)=0°. 


Figure  2.  Seven  dipole  linear  array. of  equal  x/2  spacing 
having  orientation  ((>i=O,0’=Tr/4. 
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3.  If  we  wish  to  obtain  an  electric  field  with  a Chebyshev 
expression  as 

= (cosecos;ti-sino)T7[cos(l  .2ncoS(f'Sino)]o  + 

- sin<f.  T7[cos(1 .157icos(t>sin0)](|) 


where  T7  is  the  Chebyshev  polynomial  of  seventh  degree  by  using  an  array 
with  N=b  dipoles  in  A/2  equal  distance  and  direction  as  are  shown  in 
Figure  4,  we  found 

V]  = 30.6/0°.  V2  = 4.645/10°.  V3  = 1.843/^  and  V4  = 1.0/0°. 


Figure  4.  Linear  array  with  8 non-parallel  wire  antennas. 


In  Figures  5a, b we  can  see  the  normalized  E°(()),tt/2)  and  E^(i,n/2) 
fields. 
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k 


B.  Optimization 


I 


In  Figure  6 we  can  see  the  maximum  gain, of  a linear  uniform  array 
of  nonparallel  wire  antennas  with  and  oi=2i(i-1)/N  for  N=4  to  10 
versus  the  interelement  distance  in  the  broadside  direction.  It  is 
interesting  to  say  that  in  this  case  there  are  critical  equaldistances 
in  which  we  have  maxima.  So  in  .7a  equal  distance  we  have  the  maximum 
gain  for  all  arrays. 

In  Figure  7 is  plotted  the  maximum  directive  gain  of  a circular 
array  with  tangential  wire  antennas  in  the  directions  n''=n/4, 

(i-l)/N  + r/2  for  N=4  to  10  versus  the  diameter  of  the  array.  In  this 
case  the  directive  gain  increases  as  the  diameter  decreases. 

Another  one  example  is  the  design  of. a 6-element  linear  array  in 
equal  distance  .6a  and  directions  i>’=0,  o'’=ti(  i-1  )/lC  to  provide  maximum 
gain  in  the  broadside  direction  subject  to  the  contraint  that  for  both 
E'  and  E*’  are  required  nulls  in  the  direction  (^=50°,g=90°).  We  found 
maximum  gain  0=11.559  and  the  feed  voltages 


V]  = .0887-124°. 6. 

Vp  = .2467-175°. 5, 

Vt  = .65767-250° 

V4  = l/m. 

V5  = .285740°, 

Vg  = .42597-34°. 6 

In  Figure  8 we  can  see  the  normalized  and  E'^(<(,i/2) 

fields. 

As  a last  example  the  same  6-element  array  has  maximum  gain  0=6.767 

subject  to  the  constraints  that  (i)  nulls  are  required  to  direction 

{ i'=60°,i'=90°)  and  ((^=50°  ,0=90°)  for  the  E®  and  El’  fields  correspondingly 
and  (ii)  E'^(40",90")  = l/5  E*^(90°,90”) . In  Figure  9 we  can  see  the  normal- 
ized E*M.i),7i/2)  and  E'l>(it',7T/2)  fields.  The  corresponding  feed  voltages 
in  this  case  in  absolute  values  as  they  were  taken  from  th'’  computer 
output  are: 

V]  = 57.578  + j38.129,  V2  = 81.02  + j29.205,  V3  = 91.412  + j52.613 

V4  = 100.834  + j30.476,  V5  = 61.564  + j6.619,  = 44.851  + j6.3057 
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MAX.  DIRECTIVE  GAIN 


I 


Figure  7.  riaximum  directive  gain  of  a circular  array  with  tangential 
wire  antennas  in  the  direction  '>i  = T/4,({''  = ::/2+2(i-l  )/N  for 
N=4  to  10  elements  versus  the  diameter  of  the  array. 


21 


VI.  CONCLUSION 

As  a result  of  the  foregoing  discussion,  we  have  seen  that  the 
orthogonal  method  can  be  applied  to  an  arbitrary  array  of  wire  antennas, 
and  that  it  can  solve  synthesis  and  optimization  problems  with  good 
accuracy,  apart  from  rendering  it  easily  solvable  with  the  aid  of  a 
computer.  The  method  can  be  applied  for  any  case  of  parallel  or  non- 
parallel wire  antennas,  while  the  parallel  wire  array  is  a special 
case  of  an  arbitrary  array. 

For  optimization  problems  the  orthogonal  method  has  the  advantage 
of  using  easily  applied  formulas  without  the  necessity  of  inverting 
matrices.  But  in  the  case  of  optimization  under  the  constraint  that 
one  index  has  a given  value,  the  matrix  method  must  be  used.  That  is 
because  in  nonlinear  constraints  the  formulas  given  by  the  orthogonal 
method  will  be  more  complex.  At  least  this  problem  needs  to  be  studied 
further.  In  the  other  cases  (linear  constraints)  a comparison  between 
the  Matrix  and  the  Orthogonal  method  reveals  that  the  second  one  is 
faster  and  needs  less  computer  time  than  the  first  one.  So,  both  methods 
have  advantages  and  disadvantages  which  a designer  must  keep  in  mind. 


APPENDIX  I 


From  Equation  (17)  we  have  found  that 


‘ N1  N1 

Kin  = ^ I Y_.Y*.  S„_ 

n=l  m=l 


(I-l) 


From  the  expression  (16)  we  can  find  that  as  a function  of  the 
array  geometry,  so 

f tr  rZr 


■rr  fcr 

Smn  = + E E*  ]sino  dedd-  . 

in  Jn  10  me  nj)  m4> 


By  using  Equations  (11)  and  (12)  and  integrating  we  take 


Snm  = sino"’sini"(/\„P+Bn,p)  + sinof^coso"  C|;^  + 
+ si 00*^00$ g"*  C"p  + cos0\os0'^  Gp,p 


(1-2) 


(1-3) 


where 


J ( Kt"  ) 

2.cos[2(  i™-*„)  * A*"]  [f 


(1-4) 


r— Jl/o(Krmn)  rr  33/2(Krp^n) 

B„„  = 2.cos(AA  * 


I (XR^n) 


2 

(K-'mn)®'' 


'^5/2^**''"mn) 


(1-5) 


(1-6) 
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(1-7) 


t 


when  m=n 


f 


*^3/2 


3/2 


- J ,(Kr  ) 
- ^ "" 


( ) 


5/2 


r _ ^ 

^mn  3 


By  knowing  the  S^p  the  K^j  is  known  from  Equation  (I-l). 
In  Equations  (4)  - (7)  the  involving  parameters  are: 


'■mn  = J<V*n)^*(V''r)^*!VZ,.)^ 

^mn  “ 


Li 
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APPr.\'DIX  II 

INTLf'.RAriONS  INVOLVED  ON  SYNTHESIS 


For  the  case  of  the  synthesis  as  we  can  see  from  Equation  (25)  one 
Jouble  integration  is  applied.  The  electric  field  which  is  desired  may 
or  may  not  be  one  function  with  a closed  form  mathematical  expression. 

In  all  cases  the  description  of  the  electric  field  may  be  done  Dy 
sample  points  which  we  will  show  as  follows. 


The  inner  product  of  Equation  (25)  is  an  expression  of  the  form 


Li  = 


rCT] 

j ^^(0  ^ E'^(d,e)tj/t*(';,,o)]sine  d^d'''  (II-l) 


0 JO 


A suitable  formula  for  this  double  integration,  as  well  as  for  a single 
one,  was  found  by  Sahalos  [29-30]  by  the  help  of  Chebyshev  polynomials. 


If  we  have  the  integration 
fb  rd 

j f(i;).0)d())do 

a J c 


(II-2) 


vie  can  modify  the  limits  of  integration  as 


d = i [(a+b)  + (b-a)q] 


(11-3) 


0 = [(d+c)  + (d-c)p] 


and  have 

J = Ib-ajjd-£j  I f f(p,q)-dp-dq  . (11-4) 

^ J-1  J-1 

With  sample  points  the  Chebyshev  polynomials[30]  gives  an  approxima- 
tion to  Equation  (11-4)  of  the  form 
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ii 

h 


J 


(b-a)(d-c) 

4 


ii 

.2 


n/2  n n/2 

I"  I I" 


k=0  s=0  {,=0 


^-,2kTis  

cos cos 

n 


2mTit. 

n 


■(4k^-l)(4)l^-l) 


fcos— ,cos— ') 

(II-5) 

V n n y 

{z"  means  that  the  first  and  last  coefficients  of  the  sutr,  are  multiplied 
by  1/2). 

Expression  (I1--5)  is  written  as 


J i ^ (b-a)(b-c)  I"  5;"  A(s,e)-f(cos^,cos^^  (II-6) 

n*^  s=0  ii=0  ^ ^ ^ 


the  A(s,e)  is  the  weight  of  the  point  with  q = cos(Trs/n)  and  p = cos(nS./n). 

By  using  8^  points  for  example  we  can  make  a table  with  A(s,5,) 
which  by  putting  to  the  memory  of  the  computer  can  have  more  faster 
the  integration.  In  the  following  table  we  can  see  the  A(s,e)  for  N=8. 


Some 

examples 

to  show  the 

accuracy  i 

of  the  method  follow. 

1.  The  i 

ntegral 

J = 

,5.2 

,3.2 

dydx 

4 

■ 2 

xy 

has 

an  accurate  val 

ue  equal  to 

0.123312. 

Application  of 

the  formula 

(55) 

gave 

the 

following  table: 

n 

6 

0.123 

312 

156  4 

8 

0.123 

312 

156  484 

10 

0.123 

312 

156  484  11 

as  we  can 

see 

for  n 

=8  we  have  an  approximation  with  12  s 

ignificant 

figures. 
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2.  Ttie  intearal 


f"  1 , , 

J = exp  j2n(-  sinucos^  + i sin  '^fni  + cosu  sin-;  df'dt 

Jo  Jo  2 6 

has  an  accurate  value  of  -2.247888336+jO.  Application  of  the  formula 
( 1 1-5)  gives 


n 

...  1 

J 

6 

-2.2478  +j  0. 

8 

-2.247  888  +j  0. 

10 

-2.247  888  336  1 +j0. 

3.  By  using  Simpson's  rule  the  integral 


I ? .9 

J = sin  X sin^y  dxdy 

n Jl 

gives  seven  significant  figures  for  n=32.  ‘’’he  exact  value  of  the 

integral  is  1.68267136378.  Our  fomiula  gives: 


n 

J 

6 

1.682  681 

8 

1.682  671  36 

10 

1.682  671  363  781 

Thus  convinced  that  the  approximate  integration  of  the  function 
presents  a high  degree  of  accuracy,  we  went  on  to  the  synthesis  of  the 
antenna. 
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